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. The interaction of amphiphilic polymers with small colloids, capable to reversibly 

O 

O , stick onto the chains, is studied. Adhesive small colloids in solution are able to dy- 



namically bind two polymer segments. This association leads to topological changes 
in the polymer network configurations, such as looping and cross-linking, although 
the reversible adhesion permits the colloid to slide along the chain backbone. Previ- 
ous analyses only consider static topologies in the chain network. We show that the 
sliding degree of freedom ensures the dominance of small loops, over other structures, 
giving rise to a new perspective in the analysis of the problem. The results are ap- 
plied to the analysis of the equilibrium between colloidal particles and star polymers, 
as well as to block copolymer micelles. The results are relevant for the reversible 
adsorption of silica particles onto hydrophilic polymers, used in the process of for- 
mation of mesoporous materials of the type SBA or MCM, cross-linked cyclodextrin 
molecules threading on the polymers and forming the structures known as polyrotax- 
anes. Adhesion of colloids on the corona of the latter induce micellization and growth 
of larger micelles as the number of colloids increases, in agreement with experimental 
data. 
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I. INTRODUCTION 



The interaction between polymers and small solute molecules, such as surfactants or 
colloidal particles, attracts a great interest due to wide industrial applications and their 
biological significance^. For example, polymers are used to control the stability of colloidal 
suspensions^, block copolymer micelles are employed for targeted delivery of small colloids^, 
but essentially the interaction of biopolymers with proteins or membrane phospholipids is of 
fundamental importance^. In addition, the interactions of hydrophilic polymers and brushes 
with large colloidal particles^- or micelles^ has been a subject of intensive research. 

The presence of hydrophilic polymers induces the formation of polymer-surfactant ag- 
gregates in the form of micelles covered by polymer chains^. This cooperative association 
usually happens below the CMC of the surfactants. Large colloidal particles form aggregates 
with polymers where hydrophilic polymer chains are wrapped around colloidal particles^. 
The polymer chains can link different colloids between each other, thus resulting in the clus- 
tering of colloids^ or gel formation^. In turn, several polymer chains can adsorb on the same 
colloidal particle forming the hydrophilic corona and stabilizing it. A notable feature of such 
aggregates is that if adsorption is reversible, colloidal particles can effectively slide along the 
chain. In the case of many long chains adsorbed on the same colloidal particle, this sliding 
degree of freedom results in the formation of star-shaped aggregates with annealed number 
of arms^. 

However, the interaction of polymers with small colloidal particles has received relatively 
little attention. In this paper we precisely focus our interest on interactions of hydrophilic 
polymers with small colloidal particles that can reversibly stick to hydrophilic parts of the 
polymer backbone. In contrast to large colloids, the polymers do not adsorb onto the par- 
ticles but only adhere onto a few active sites, i.e., the small colloids induce topological 
changes in the polymer configurations, promoting looping, cross-linking, and interconnec- 
tions between the chains. In this context, the following situations, which are relevant in 
practice, can be mentioned: (i) The driving force for the formation of mesoporous silica- 
based materials type SBA or MCM— is the strong interaction of hydrophilic polymers with 
silica particles, e.g. TEOS, and the subsequent formation of self-assembled structures. The 
possibility of tuning and controlling the design of the resulting structures is a challenging 
task. Silicon alkoxides in water form Si-OH compounds, but also different pairs or oligomers 
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of silicon, namely Si-O-Si^. Silica particles form hydrogen bonds with hydrophilic poly- 
mer's backbone^ . Since the hydrogen bonds are relatively weak, the silica particles can 
effectively slide along the chain by breaking and re-forming new hydrogen bonds. Thus, 
we can consider a silica particle as being either a sliding link between different chains or a 
sliding loop on the same chain. Addition of silica particles to the solution of block copolymer 
micelles leads to the adsorption of silica particles on the coronas of micelles^. The practical 
example of such systems is an aqueous solution of TEOS with common triblock copolymers 
Pluronic P123 (E0 2 o-P0 70 -E0 2 o) or Pluronic F127 (EOi 6-P0 7 o-EOio 6 ). The micelles and 
self-assembled ordered structures of these polymers serve as precursors for the pores in the 
mesoporous materials. The diameter of hydrophobic cores of micelles determines the size of 
the pores, while the structure of the polymeric corona is responsible for the microporosity 
of the material^. Therefore, the control of the adhesion of the silica particles into the coro- 
nas can lead to the control of the core of the micelle and, eventually to the control of the 
properties of the final porous material, (ii) A cone shaped molecules with a hydrophobic 
cavity known as cyclodextrins^— can spontaneously thread on linear hydrophilic polymers. 
Polymer chains dressed with necklaces of cyclodextrins bear the special name of polyrotax- 
anes. Since two molecules of cyclodextrins can be covalently bound together, the resulting 
molecule will represent two linked rings. Such molecules can either link two chains or form 
a loop on the same chain, providing a sliding degree of freedom for the chain to move inside 
the rings. In previous works the effect of these sliding links for the polymers grafted to a 
surface^ and the micelles with sliding coronas^ has been studied. 

In this paper we introduce a formalism suitable for a qualitative description of these 
mentioned processes, in which flexible polymers interact with small colloidal particles or 
other agents, that induce binding with the polymers. A complete statistical mechanical 
description of these processes is of extreme complexity due to the changing topology of the 
system due to the adhesion of the colloids in different chain backbones and, therefore, we will 
here focus on the qualitative trends that a scaling theory can reveal. Formally, the complex 
topology of the polymeric network, produced by the presence of the binding agents as shown 
in Figure (TJ permits us to use known results of polymeric networks of fixed topology, for 
which the relevant scaling exponents are known^S. The challenging objective of our work is 
thus to derive the partition function for a system where the vertices of the polymeric network 
can slide along chain backbones and, therefore, the topology of the system is not fixed. This 
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FIG. 1. Schematic representation of a block copolymer micelle with small colloids in the corona. 
This sketch shows possible topologies addressed in this work. 

is precisely the case of the physical situations that we have described so far, and that cannot 
be trivially addressed by a simple mapping of the fixed topology results, as done in other 
context^ - — . An important consequence of this fact is that, in the integration of this degree of 
freedom, the system passes through regions with different topology, characterized by different 
sets of scaling exponents. In the next section we present what, to the best of our knowledge, 
is the first derivation of this effect. Therefore, in the following we formulate our scaling model 
and address the question of the calculation of the singular non-trivial contribution to the 
entropy of a system of hydrophilic polymers with small colloidal particles forming different 
mobile polymer architectures. With the construction of this singular part of the partition 
functions, according to the scaling theory, we find the most favorable conformations of the 
system and, as a consequence, the distribution of colloids along the chains, for instance. 

Within this formalism we can calculate the regimes of loading of silica particles in lin- 
ear polymers, star polymers, and in micelles of block copolymers, and permits us to give 
a theoretical justification of several experimentally observed facts. We believe that our 
approach can be useful in biophysical applications in systems like DNA chains interact- 
ing with proteins related to gene regulation models, where equilibrium analyses have been 
carried out^, and also the effect of looping on the transport of proteins along the chain 
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has been analyzed^ 5 .^. Related analyses can be carried out to study the problem of DNA 
denaturation2I and, effectively, have been done to study the equilibrium shapes of knots in 
chains^i2£. 

The paper is organized as follows. In section [III we introduce the model of the system 
and perform the scaling analysis. Section III II is devoted to the study of the equilibrium 
between polymeric systems and colloidal particles, in particular, colloids with star polymers 
and colloids with block copolymer micelles. Finally, section [V] is devoted to the conclusions 
that we have drawn from this work. 

II. MODEL 

In this section we propose the scaling form of the partition function for systems, poly- 
meric networks, of different topology. Such a partition function contains a trivial energetic 
contribution arisen from the short-range interaction between the binding agents and poly- 
mer segments. However, the entropic contribution strongly depends on the topology of 
the polymeric network resulting from that interaction, together with the topological effects 
that are related to the topology of the polymer itself, as in the case of star polymers. To 
that purpose, we refer to the results obtained by Duplantier— for the partition function 
of polymer networks of arbitrary topology. Such a partition function has a superscaling 
form Z ~ s N N' yc ~ l g(ni/N,n2/N, ...,rik/N), where N is the total length (total number of 
monomers) of a network, s is a non-universal geometrical constant, m, nk are the lengths 
of different polymeric threads between two crosslinks, or between a crosslink and a dangling 
end, such that n j = ^5 an d 9 is an unknown function. 7£ is the universal exponent which 
does not depend on the type of interactions and is determined only by the topology of the 
system. It is given by^ 



where qi oop is the number of independent loops, pk is the number of vertices with k legs, d is 
the dimension of the space, v is the Flory exponent, associated to the radius of gyration of a 
self-avoiding walk, is the exponent corresponding to a vertice with k -legs. The universal 
exponents v and <Jk are known exactly for d = 2 and d > 4 and numerically for d = 3. In the 
following we use the values of obtained from the simulation results for critical exponents 
of star polymers^. Although a general form of the function g(ni/N,ri2/N, ...,nk/N) is 




(1) 



fc>i 



5 



not analytically known, the scaling behavior of the system can be derived from the scaling 
form of g when some of its arguments go either to or to 1. This scaling form gives 
the correct asymptotes in these limiting cases. Due to the fact that our interest lies only 
in the derivation of the scaling behavior of the overall system, we will avoid the explicit 
construction of the crossover functions between different regimes, characterized by different 
topologies, by basically assuming that the crossover function is of the order of a constant 
and matching these constants at the crossover region between different topologies, as shown 
shortly. This method allows us to derive an approximate description over the whole range 
of topological configurations of the system. From a physical point of view N as well as 
series of ni,ri2, ■ ■ ■ , representing the size of the branches of the polymer network represent 
the number of Kuhn segments of each polymeric thread. Furthermore, the validity of the 
scaling analysis is restricted to a system with N being very large. In particular, Monte 
Carlo simulations on isolated polymer s^ 1 ^ indicate that the scaling limit is reached for 
sizes exceeding hundreds of Kuhn segments, and therefore our analysis will be restricted to 
polymeric systems of at least this size. With regard to the colloidal particle, we will consider 
that its size is slightly larger than the Kuhn segment but much smaller than the overall size 
of the polymeric network. 

We base our model on the following assumptions regarding the polymer-small colloids 
interactions: (i) one sliding sticker can bound several polymer units either on the same 
chain or different chains; (ii) the links are reversible and, hence, the reversibility allows for 
an effective sliding of the sticker along the chain; (iii) we assume steric repulsion between 
small colloids. Then the colloids and polymers bind together they define a network whose 
topology can change due to the sliding degree of freedom of the colloid. It is important 
to realize, however, that the sliding can physically take place by reversibly breaking and 
reforming bonds between the polymer and the colloid, regardless whether the colloid com- 
pletely unbinds from the chain or it slides from site to site, since both are dynamic processes 
whose effect is accounted for by the equilibrium statistical mechanical treatment. Hence the 
challenging objective of this work is to determine the partition function of a system with 
fixed number of polymers and binding agents, taking into account the degree of freedom of 
the motion of the colloid along the chain backbone and, at the same time, correctly account- 
ing for the topological changes in the network originated by this degree of freedom. Up to 
the best of our knowledge, this task has not been carried out up to date. 
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TABLE I. Numerical values of the critical exponent v and star exponents a, obtained by interpo- 
lation of Monte Carlo results for star polymers^ in d = 3 dimensions. 
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-3.709 
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-6.084 


-6.975 


-7.916 


-8.898 


-9.926 


-10.991 


-12.093 



In what follows we will apply the formalism to systems of increasing complexity, empha- 
sizing the details of the calculations in the simpler systems and making the extension to 
more complex systems. 

A. Stickers in linear chains 

Let us first examine the case of an isolated chain with a single sliding sticker binding two 
Kuhn segments, as shown in Figure |2j In this figure, one can identify one crosslink and two 
chain ends, which divide the chain into three threads. If the chain has a total length N, the 
first thread is a tail of size I, the second is a loop of size n, and then another tail of length 
N-n-l. 




FIG. 2. A single loop of n monomers created by a sticker at a distance I from one of the ends. 

The partition function of the structure sketched in Figure [21 has a general form Z ~ 
s N N' rc ~ 1 g(l/N,n/N, (N — l — n)/N), where the exponent is given by eq. (P), with qi oop = 1, 
with two vertices of one leg a\, and one vertex of four legs 04, i.e. 7£ — 1 = 2a 1 -Va^ — dv (see 
Tabled for numerical values). Without the explicit knowledge of the crossover function g, 
the exact partition function for the system of Figure with the sliding degree of freedom, 
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cannot be obtained. However the scaling behavior can be derived from the singular form of 
g in the appropriate limit, obtained by making some of its arguments going to or to 1, 
where g is expected to behave either as a power law or as a constant. With this in mind, we 
propose a power law form of the partition function that will permit us to properly interpolate 
between the singular behavior of the crossover function g in the different limits, i.e. 

Z ~ n x l v (N — n — l) z (2) 

where the exponents x, y, z can be fixed from the known limits. In this simple case, one of 
these limits is, for example, l/N — > 0, n/N ~ 1, and (N — n — l)/N ~ 1. Our procedure will 
become apparent in what follows. 

The partition function of the system requires the integration over the sliding degree of 
freedom. The final result has a general structure of the form 

Z(N, m) ~ e~ m£/kT s N Z(N, m) (3) 

where m is the number of stuck colloids, and the factor e - m£ / fcT takes into account the 
interaction energy in the sticking process. The entropic part contains the trivial factor 
^depending only on the number of monomers contained in the polymeric network. Z takes 
into account the nontrivial contribution to the configurational part of the partition function 
which depends on the topology of the polymeric network, but also on that induced by the 
presence of m colloids adhered to the polymer. From now on, we will concentrate on Z, and 
omit the tilde, for the ease of notation, where confusion could not occur. 

For the case of one polymer chain and one colloid, the nontrivial part of the partition 
function takes the form 

pN/2 pN-21 

Z(N,m = l)~2 dl dng{l/N,n/N,{N -l-n)/N) (4) 

Jo Jo 

The integration limits take into account that the system is invariant under the permutation of 
the two tails and therefore we can restrain the integration domain by the line I = (N — n)/2, 
which implies that the tail of size I is the smallest tail. We have to further split the integration 
domain into separate regions of different topologies due to the fact that each topology 
requires a set of exponents x, y, z in equation ()2] ) to consistently interpolate between the 
appropriate limits. We can argue that these regions are characterized by I) small loop, 

8 



n < I < N—l—n, II) intermediate loop, I < n < N—l—n; and III) big loop, / < N—l—n < n. 
Due to the fact that we take I as the smallest tail without loss of generality. We when have 



i-N/3 pi pN/2 rN-21 

Z(N, m = 1) ~ / dl dnn x l v (N - n - l) z + dl dnn x l y (N-n- l) 1 

Jo JO J N/3 Jo 

pN/3 rn 

+ d / dn dln x 'l y ' (N-n- l) z> 
Jo Jo 

r-N/2 pN-2n 

+C X / dn dln x 'l y ' (N-n- l) z ' 

Jn/3 Jo 

N/2 AN-n)/2 



+ C 2 dn dln x l y (N-n- l) z 

Jn/3 Jn~2u 

pN p(N-n)/2 

+C 2 / dn dln x "l y " (N-n- l) z " (5) 

Jn/2 Jo 

where the constants C\ and C 2 measure the relative weight between the functional form 
taken for each of the contributions. In what follows we will determine the three sets of 
exponents and evaluate the contributions of each term to the partition function. 

II) Intermediate loop, I < n. In this region, only three limits are possible. In the first 
limit we consider that all legs and the loop are of a similar size and of order N. Therefore, 
according to eq. (ED, Z ~ N x+y+z . Then, making use of the general expression for the 
topological exponents given in equation ([T]), one can write 



x' + y' + z' = 2<7i + 0-4 - du (6) 

corresponding to a topology of two free ends at the extremes of the two tails, plus a vertex 
of four legs where the loop closes in the middle of the polymer. In the second limit, if the 
small tail vanishes, I — >■ 1 (l/N — >■ 0) with n ~ N/2 and (N - n - 1) ~ N/2, the topology of 
the system evolves towards a long tail with the loop at the end, characterized by a topology 
with the presence of one free end at the extreme of the remaining long tail, and a vertex of 
three legs for the loop closed at the other end of the tail. Then Z ~ n x 'l y ' (N — n) z ' . Since 
this expression has to be compatible with the known singular behavior of the system with 
the mentioned topology ([!]), we further have 



x' + z' = <j\ + 03 — du 



(7) 



Finally, in the third limit, if both, the small tail and the loop, vanish, I < n — > 1, we get a 
free chain topology and thus 



z' = 2cj x (8) 

These three equations permit us to find the numerical value of the desired exponents for the 
interpolating function that will be used in the integration to derive the singular behavior of 
the partition function. We obtain x' = 03 — o\ — du, y' = 0\ + 04 — 03, and z' = lay. 
We can now write the partition function of this domain in the form of nested integrals 

Z n ~ / dnn a3 ~ ai - du / dir +<T4 - a3 (N -n- l) 2ai + 
Jo Jo 

pN/2 rN-2n 

+ / dnn az ~ ai ~ du \ dll ai+a4 ~ a3 (N-n- lf ai (9) 
Jn/3 Jo 

In both terms, since 0\ + — a 3 —0.2 > —1 and 2cxi = 0.16 > 1—, the inner integral 
converges. However, the inner integral in the first contribution behaves as n <J1+ai ^ az+1 
for small n. Therefore the outer integral is diverging due to lower bound, leading to a 
contribution of the order A <T4-dl ' +2 iV 2<Tl w A CT4_ ^ +2 A r °' 16 , where A is the lower cutoff for 
the size of the loop, of the order of one monomer. The second contribution is convergent, 
yielding N a ^ du+2ai+2 dtt a ^- dv J^~ 2t dss ai+a4 - a3 (1 - t - s) 2ai = o.31A^-^ +2cti+2 w 
0.31iV" 098 . Therefore, in the limit N — > 00 the partition function is dominated by the first 
term 



Z H ~ A a4 ~ d ' /+2 N 2 ' 71 (10) 

This result indicates that big loops are not entropically favorable and that the system will 
tend to reduce the size of the loop to the minimum. Moreover, since in this limit the small 
tail is, by construction, smaller than the loop, the most probable conformation is hence a 
minimal loop at the end of a long tail. 

I) To evaluate Zi, one has to consider that the loop is smaller than the shortest tail, 
n < I. When the sizes of the loop and the tails are comparable, we have the same scaling as 
in the II-part. Hence, in the latter case 



x + y + z = 1a\ + <7 4 — du 



(11) 
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Notice that when the loop and the tails are comparable and of the order of N, the integrands 
of the II- and I- part should be comparable. Since the functional form is also the same for 
both cases, we have that the constant C\ ~ 1. 

As before, to fix the exponents we identify the singular limits compatible with the con- 
dition n < I. Taking n — > 1, the free chain limit is recovered and hence 



y + z = 2a 1 (12) 

Furthermore, when the small tail vanishes, this implies a simultaneous vanishing of both the 
loop and the tail, n < I — > 1. One gets 



z = 1o x (13) 
Finally, one also obtains x = 04 — dis, y = 0. The partition function reads 



N/3 fl i-N/2 pN-2l 



dl dnn^- du (N -n-iy ai + dl dnn ai (N — n — I) 1 (14) 

JO J N/3 JO 

Since again the integrand of the inner integrals diverge the integral is dominated by n — » A. 
The outer integrals have a regular integrand and the whole term scales as jV 2<71+1 . Then the 
final expression for the partition function of this part is 



Z/ ~ A CT4 "^ +1 iV 2CT1+1 (15) 

This partition function indicates that under the condition n < I the most favorable con- 
formation is the existence of a minimal loop traveling along the chain. Notice that this 
contribution has an extra N factor as compared to the case //, that is Zj/Zu ~ N/A — > 00, 
for very long chains, indicating that this contribution will be dominant. Similar conclusion 
is drawn for the case knots in chains^. 

Ill) The third contribution in (jSJ) is obtained as before. For this case the relevant limits 
are: both tails and the loop are of the same order, one tail vanishes, I — > 1, both tails tend 
to vanish I < N — I — n — > 1. With these limits we find x" = —dv , y" = <j\ + 04 — 0-3, and 
z" = G\ + 03. Furthermore, since Zm should be comparable with Zi and Zu when the two 
tails are of the same order we conclude that C 2 is of order 1. 
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N p(N-n)/2 

dn / dln- du l ai+ai - (73 (N - n - l) ai+a3 + 

N/2 JO 
N/2 r (N-n)/2 

dujai+CT4— (73 / a j ;\o-i+<T3 



'111 



dn / dln- dv r +aA ~ aA (N -n-lf 1+(7Z (16) 

7V/3 JN-2n 

Both integrands have integrable divergences and, therefore, the N dependence of Zm can 
be trivially obtained as before N a *- du+2ai+2 [/* dtt~ dv J (1 "* )/2 dss<™- a * (1 - t - s) ai+,T3 



44jy°"4-a!^+2(Ti+2 



+ jljl dtt- d » f^_~J /2 ds8 n+ "*-°* (1-t- Sp +U3 

We can then conclude that the dominant contribution to the partition function of a 
sticker in a free chain is given by Z ~ Zj ~ ^0-4-^+1^20-1+1^ s j nce Zu/Zj ~ N~ l and 
Z Jir /Z 7 ~ ^04-^+1^04-^+1 _ TV -126 . 

Therefore the overall iV-dependence of the partition function is that of a bare chain with 
an extra power of N that takes into account the freedom of the location of the small loop 
along the chain. The loop cut-off size A could be linked to the local stiffness of the chain 
(persistence length). The optimal configuration is a small loop freely sliding along a linear 
chain. This is a general result and, as we describe later, is applicable to the loops sliding in 
the corona of a micelle or the arm of a star polymer. The position of the loop I does not 
enter the final expression and the distribution of positions of the small loop along the chain 
and hence the probability of finding the loop (the sticker) at a position I along the chain 
backbone P(l) is uniform up to a distance A from the chain ends. 



N-a-l-n r n 2 




FIG. 3. Ailed loop created by two stickers on the same chain. 



To end this section, let us consider the case of m stickers in one linear chain to the light of 
the results that we have just derived. As we have seen, the more important contribution to 

12 



the partition function is the degree of freedom of each small loop along the chain backbone. 
This contribution is lost if the loop becomes large as well as in conformations in which loop- 
like structures are linked by a sticker. For example, in ailed loops (Figure [3]) the translational 
entropy gives a factor N while, if the loops were independent, the entropy would be of order 
N 2 (see Appendix [A]). Then the partition function for a system of m stickers in a free 
chain can be straightforwardly generalized considering the independent contributions of the 
individual loops, provided that mA <C N. Thus, from the expression (TT5]) one obtains 

Z ~ /\rn{a 4 ~dv+l) N m N 2a x 

Notice that the factor N 2ai stands for the entropy of the supporting chain. The result 
given in eq. (JTTjl is only valid in the case that the sticker-sticker interactions are negligible 
and therefore, when the chain is not saturated. 

B. Stickers in star polymers 

In the following, let us assume that the star polymer has p arms, each of length N. 
We will proceed by analyzing the dominant diagrams starting with one colloid and then 
generalizing the result to an arbitrary number m of stickers adhered to the arms of the 
polymer as before. 

1. One loop on one arm 

The simplest case is that of one loop in one of the arms of the star polymer, as shown 
in Figure HI The main difference with the respect to previous analysis is here the presence 
of the topological constraints at the vertex where the p arms meet, giving strong excluded 
volume interactions in the vicinity. This is reflected in the partition function of the system 
by the presence of highly negative exponent a p corresponding to strong excluded volume 
effects in the center. 

Denoting the size of the loop by n and the distance of the loop from the center by Z, 
we write the partition function of the system star polymer - colloid under the interpolating 
form already used Z ~ n x l y (N — n — l) z , where x, y and z are exponents to be determined by 
matching the appropriate limits. Furthermore, as before, the set of exponents is not unique 
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FIG. 4. Sliding loop of size n on the arm of a micelle at the distance I from the center. 

but depends on the relative sizes of the three entities in which we have divided the chain, 
i.e. loop, tail, and distance to the vertex. We demand that these different expressions for 
the interpolating function continuously crossover at the limits of their validity. 

Following the same analysis of the previous section, the interpolating function have to be 
split into six different sets of exponents due to lack of symmetry of this case. However, we 
know that the partition function is going to be dominated by a small loop traveling along 
the arm. To simplify our analysis we will concentrate on the scaling form of this regime, 
defined by a set of exponents in the region n < I < N — n — I, and a different set in the 
region n < N — n — I < I. After analyzing all the corresponding limits, one has 



pN/3 A r-N/2 [-N-21 

Z~ dl I dnn ai - du (N -n-l) ap+pai + / dl dnn Ui ~ du {N - n - l) a v +pui 

Jo JO J N/3 Jo 

i-N/2 p(N-l)/2 i-N f{N-l)/2 

+ / dl j dnn a4 - du l^ +pai + / dl dnn a4 - du F p+pai (18) 

J N/3 J N-21 J N/2 Jo 

As expected, the inner integrals diverge for this small loop case when n — > A, which 
reflects the aforementioned dominance of small loops with respect to large loops. After 
integration we get that the dominant contribution scales as 
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Z ~ /\°4-du+l N a p +p<n N ^g^ 

The term N r7p+pai is the partition function of a star polymer of p arms, while the additional 
N factor stands for the translational entropy of the sliding small loop. The other topologies 
that we have not explicitly analyzed imply either a loop of order N or small loops that are 
interacting with the vertices. In all these cases the translational entropy of the loop is lost 
and these contributions are therefore subdominant. The integrand in the partition function 
(j!8p . upon integration over n, allows us to derive the distribution of these small loops along 
the chain, giving 

P(l) ~ ( K ' ' ' (20) 

l a *+P' T \N/2 < I < N 

which indicates a slight repulsion of the small loop from the center of the star. 



2. Several loops on the same arm 

The partition function of several loops on the same arm can be obtained by a straight- 
forward generalization of the previous result, as we have done for the single chain case. We 
obtain 



which basically contains the bare entropy of the star polymer together with N m factor due 
to the translational entropy of m loops. We have implicitly considered that the number of 
monomers in the z^-branch, rii is of order N. This result is however limited to the case 
mA N, due to the fact that no interactions between loops has been considered. 

This result can be further generalized to the case where the m stickers are distributed 
among p branches. Let us consider a star polymer with p -arms of rii, n 2 , . . . , n p monomers 
each, with Yui=i n i = ^> where colloidal particles can adhere from a given bulk solution. 
The complete partition function of m colloids attached to this polymer is given by 

Z(m) OC £ e -ms/ k T s N Am ^ + l )N m ^ ^ ^ 

mi +7712 H \-m p =m 

(22) 
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where use has been made of eq. (l2Tj) to introduce the dominant contribution to the entropy of 
the system, and Y^i=i m i = m - Branches are assumed to be discernable due to polydispersity 
considerations. However, if the size of the branches n; are comparable and rij ~ N n\/N ~ 
712/ N ~ . . . ~ n p /N we can approximately write 

Z ( m ) oc ( m +P-y i e -me/kT s N A m(« 4 - d „ + l) ( N\ ™ 

m\(p — 1)1 \P J 

The combinatorial factor stands for all combinations of mi, wi2, . . . , m p such that its sum 
gives m. 

3. m stickers looping in many arms 




FIG. 5. A loop between two arms of a micelle created by one sticker in the center of the micelle. 

The exact consideration of the general case is a formidable task in view of the multiplicity 
of exponents and integration domains that we have shown in the simple cases analyzed 
above. However, we prove here that when two arms are connected by one sticker, there is 
a significant entropy reduction with respect to the case where the sticker form one loop in 
one arm. Therefore, the dominant configuration for the general case is that in which the m 
stickers form small traveling loops in individual arms. To prove this, let us first consider a 
loop formed by two different arms as shown in Figure O The sticker splits a first chain into a 
segment a close to the center and a free tail N — a. The second chain is split into a segment 

1(3 



b and a free tail N — b, while the loop size is n = a + b. Since two arms are symmetrical we 
sort the tails by their respective lengths, a < b this leads to N — b < N — a. The general 
form of the partition function is 

Z ~ a Xl a X2 (N -a) Zl (N -bf 2 (24) 

where the exponents x±, X2, Zi, z<i are fixed by the limits. The dominant contribution to 
the partition function is found in the limit a < b < N — b < N — a. This can be intuitively 
understood by noticing that the size of the loop is directly related to the location of the 
sticker. In this way, when the sticker is located near the chain end (a large loop) there is a 
strong entropy reduction due to the reduction of the two tails to create a big loop. Hence, 
the dominant contribution is that of a small loop near the core of the star. This argument 
also indicates that the translational entropy of the small loop discussed in the previous cases 
is here lost. Therefore the dominant contribution to the partition function reads 

pN/2 pb 

Z~ dbb**-*-**-*' \ daa a ^+ 2+a4 (N - a) Zl (N -bf 2 (25) 
Jo Jo 

where use has been made of the relevant limits to fix the exponents, with Z1 + Z2 = cr p +pcri 
as imposed by the limit a, b — > 0. Since a p — a p+ 2 + cr 4 > —1, the inner integral converges, 
giving, Z ~ j^ 2 dbb a4 ~ dl/+l N crp+I " 71 , while this integral diverges on the lower limit A, since 
<jp+2 — cr p — dv < —1 for all p. One finally obtains 

Z ~ ^o-4-^+2 jy-<7 p +p<7i (26) 

Similar arguments lead to the conclusion that a larger number of loops between two 
different arms or nested loops of several arms as shown in Figure will degenerate into 
immobile loop in the center: lower limit exponent < — 1 for any p, including p = 4, and the 
nested integral in (125]) is of order of b regardless the integrated function. 

III. EQUILIBRIUM BETWEEN COLLOIDAL STICKERS AND A 
p- ARMED STAR 

With the help of the results obtained so far we analyze here the equilibrium between a 
solution of small colloidal particles and polymeric structures. This method can be applied to 
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FIG. 6. Two examples of nested loops formed by two stickers in the arm of a micelle. 

determine the partition function of a system in which small colloids can adhere to polymers 
of complex topology. On one hand, it is useful to find how the topological constrains can 
modify the preferred location of the colloids in the system, but also how colloids can induce 
the binding of different branches. Therefore, it is not a priori evident the final equilibrium 
structure. 

In the first place, we study the case of a star polymer of p-arms in equilibrium with a 
suspension of colloids. In the second place, we will drive our attention to the case of block 
copolymer micelles, where the equilibrium can change the properties of the micelle as such. 
This latter case is of interest in the moulding of porous materials from the self-assembly of 
surfactants through silica colloids associating with the corona of copolymeric micelles, as it 
is the case for porous materials of type SBA or MCM, for example. 

Making use of eqs. ( |23i) . the grand partition function for the system then reads 



S oc J2 e" m/kT Z(m) ~ s pL {pLf^ ]T [e^ A^ d+1 L~\ ™ ^yf^W (27) 

m=0 m=0 

Stirling approximation can be used if m » 1, while p can be of order 1, thus the combi- 
natorial factor turning into m p ~ 1 / (p — 1)!. Furthermore, we have introduced the number of 
monomers per arm, L, so that N = pL. The grand-canonical potential Q is defined as usual 
f2 = — ZcTlnS. Defining Y(fi) = e 1 ^ A a4 ~ ud+1 L one finds that the sum is converging always 
that Y(ji) < 1. In particular, in the limit Y{p) — > 1 we find 

2 x (T^RF (28) 
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Therefore, colloids saturate the star polymer when the chemical potential rises the critical 
value fi c = e + kT\n (A a4 ~ ud+1 L) . In particular, we find that the average number of colloids 
that aggregate in the star is given by 

(m) Y.ZvYfjj.rm* i r(p + i) + c(-P) „ P (2q) 
{ ) YZ^^r^ 1 i-y(j*) r(p) + c(i-p) ~i-y(jm) 1 J 

where the last expression corresponds to the limit Y — y 1~ and p is taken p > 4. Therefore, 
above the critical value of the chemical potential the star saturates and m becomes of the 
order of 1/a ^> 1, where a is A/N = A/pL. The quantity am can be regarded as the 
fraction of polymer wrapping colloids, which tend to 1 at saturation. 

The analysis of the saturated state cannot be done with the same rigorous approach 
developed so far. We will assume that saturation occurs because of the limitation of the 
available polymer for the colloids in one arm due to the presence of other colloids in the 
same arm, in the spirit of van der Waals equation of state. One can then approximately 
write (L — mA/p) m instead of L m and (L — mA/p) crp+pcri m eq. (]23|) . to account for the 
limitation of the size of the polymeric threads at fixed number m of minimal loops. This 
approximation will be valid always that mA <ti N. One can finally write 

p-l / A \ m / a \ (Tp+pa! 

Z(m,p) oc e -™/kT s pLA m ^- d " +1) L m (p L)^ +pai 1 - m— 1 - m— 

(p-l)! V PLJ V pLJ 

(30) 

The partition function can then be written as 



N / A / m A \ m + CT P+P°"i 



m=0 

Let us center our attention on the value of the sum 



E^r(l-^) mT 1 (31) 



1/a 

/„_! = £ F(/i) m (1 - am) m+a * +pai m"" 1 (32) 



m=0 

l± 



in the limit a — > 0, near the transition, Y — > 1 . Then, the summand can be rewritten as 

y(/i) m (1 — am) m+C7p+Pai m P ~ 1 ~ e mlnY+(m+a p +pa 1 )ln(l-am)+(p-l)him ^gg^ 

The exponent in the second expression can be regarded as an effective Hamiltonian 

H{m) = — (m In Y + (m + a p + peri) ln(l — am) + (p — 1) In m) (34) 
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To analyze the stability of the system, let us calculate the saddle point m* by demanding 
that 

= (35) 



OH 

dm 



m* 

This gives the so-called equation of state, which reads in our case 

InY + ln(l - am*) - ^±^±^j + 1^1 = (36 ) 

1 — am* m* 

Let us assume that near the transition, InY <C 1, the fraction of wrapped polymer am* is 
small. Then, from eq. (I3T)|) we obtain an estimate of the average value of m*. To the lowest 
order in a we find In Y — 2am* — a(a p + pa±) + (p — l)/m* = 0, which gives 

m* = lnY - a{ ^> + Pai) ± ^8a(p - 1) + (hi Y - a(a p + pa,)) 2 (37) 

The negative sign gives an unphysical m* < solution. This expression gives a crossover 
between values of Y smaller than 1 and values larger than 1 for finite a. The assymptotic 
behavior is obtained by expanding the expression for small a, in fact for hiY/yfa larger 
than 1. At the crossover, |lnY| ~ a 1 / 2 . Therefore, if 1 > | \xiY\ > a 1 / 2 the solution reads 

InY p-1 2(p-l) 2 



m 



2a In Y In 3 Y 

p-1 2(p — l) 2 

lnYl ~ |ln 3 Y| 



-a + O {a In Y) for Y > 1 (38) 



m --nr^~ m 3,,i a + 0(akiY) for Y < 1 (39) 



For values of the control parameter a 1 / 2 > \ \nY\ > 0, the appropriate expansion reads for 
both branches 



m 



p_l_e_ L +_pa 1 + W f w 2 lnY\ 
2a 4 4a V a 1 ' 2 J 

This result indicates that the free energy is analytic through the transition and the saturation 

is therefore a crossover process as the chemical potential increases for finite a, due to the 

existence of this crossover regime a 1 / 2 > \ \nY\ > 0. However, the range of validity of eq. 

PU|) decreases as a — > indicating that the behavior of m* will become singular precisely in 

this limit. To further analyze the character of such singular behavior, we have to consider 

the limit InY — > with |lnY| > a 1 / 2 . The appropriate order parameter in this singular 

limit is the wrapped fraction of polymer, ip = am*. Since for ip = am*_ = for Y < 1 while 

ip = am* + ~ In Y if a — > 0. In addition, notice that am* ~ a 1 ! 2 — > in the crossover regime, 

eq. (T4*0|) . Let us further define ( = In Y 1 . We then rewrite the effective Hamiltonian as 



aH(m)-a(p-l)\na = h(ip) = - (ip( + (ip + a(a p +po"i)) ln(l - ip) + a(p - l)ln^) (41) 
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Then, the main integral eq. ( )32|) reads 

1 f 1 , , hW 



Ip-i-— I d^e~— (42) 



o 



The saddle point analysis then yields 

P-1 2 2(P~1 
— h a — a 

2 C C 



^.~i + a^ a 2 KF A3 ; + OC) forC>0 (43) 



^* ~ - a 2 2 ^ 3 ^ + O (aC) ~> for ( < (44) 

where we recall that £ > a 1 / 2 . Furthermore, the mean field stability of the solution requires 



<9# 



a^— F + 2 + a(<7 p + pa x ) + (3 + 2a(<7 p + p<7i))^* + O (^*) > (45) 
where second equality follows considering that we are near the transition and therefore 



« 2 + 3^* + O (^*) 3 > (46) 

1> 



xjj <C 1. Notice, however, that the limit a — > has to be taken before x/i — > 0, i.e. before 

we approach the transition. Therefore, Since h ; '(ip) > in all the range of values of the 

mean field solution is stable. Therefore, 

d 2 h 
dip 2 

which is always positive, indicating the stability of the mean field solution. Near the tran- 
sition, the energy is continuous and behaves as 

h(r+) ~ -j (47) 

h(ip*_) a(p - 1) sign(C) -> (48) 



The function h"(i(j*) is also continous at £ = 0, since 

h"(r+) ~ 2 + > (49) 
~ 2 > (50) 

In conclusion, for a — > at finite |£| <C 1 the system presents a crossover from a non- 
decorated state, = 0, to a decorated state, ~ £/2, in which the polymeric threads 
are increasingly covered by colloids. However, it is noteworthy that the second derivative 
of the saddle point free energy is continous in the transition, but we observe that the third 
derivative has indeed a finite jump. Therefore this behavior can be interpreted as a higher 
order phase transition in the sense of the old Ehrenfest classification, but not in the sense 
of an order - disorder, Ising-like transitions. 
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IV. MICELLES OF POLYMERS WITH STICKERS 



Colloidal stickers can aggregate in coronas of block copolymer micelles and modify the 
equilibrium structure of micelles, which is precisely the relevant problem for the formation 
of mesoporous materials through self-assembly of block copolymers. After the discussion of 
the saturation transition, we address this structural problem. With this purpose, let us first 
construct the free energy of a micelle with p block copolymers and with m adhered colloids 
to consider afterwards the global equilibrium of a system of many different- sized micelles 
with free chains and colloids in the bulk. 

Since the colloidal stickers interact only with the hydrophilic blocks forming the corona, 
only the corona contribution is affected by the presence of stickers and the partition function 
can be split into two factors, namely, the contribution due to the corona and that due to 
the core. 

Thus, the free energy of corona F corona is given by 

Fcorona(m,p) = - hxZ(m,p) (51) 

where Z(m,p) is given by (1301) and we have explicitly shown the parameter p that will be 
relevant in the following discussion. Implicit in eq. f[5"Tj) is the fact that the size of the 
hydrophilic block is much larger than that of the hydrophobic, so that the micelle can be 
regarded as effectively a star. 

This repulsive contribution of the corona, which tends to solubilize the copolymers is 
balanced by the attraction of hydrophobic units in the core. The core contribution F core (p) 
can be written as 

Fcoreip) = ^Rl{p)l (52) 

where 7 is the surface tension between the core and the solvent, being the core of size R c . 
Assuming dense packing of the monomers in the core, the radius of the core composed by p 
block copolymers can be expressed as 

1 

(53) 

where v is the volume of a monomer, L c is the length of the hydrophobic block. Finally, a 



Rc(p) 
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contribution, F con f(p) should be added due to the entropy reduction due to the aggregation 
of p chains into a micelle, yielding 



F C onf{p) = (p- l)ln 



p 

V 



(54) 



See Appendix [B] for the details. We consider that the core is compact and thus V agg ~ R^; 
hence V agg ~ vpL c . 

Hence, the free energy of a micelle of p block copolymers and m stickers in the corona is 
the sum 

^corona 

(m, p) + F core (p) + F conf (p) (55) 

where p > 1. The particular case F(0,p) is the free energy of a clean (with no colloids 
adhered) micelle of p block copolymers, while F(m, 1) is the free energy of a free chain 
decorated with m colloids. To establish the equilibrium, we have to further introduce the 
free energy of isolated colloids which read, F coUoid = cyln(^A^), where cj stands for the 
concentration of free colloids in the bulk and is the de Broglie's length of the colloid. 
Notice that we consider that the bulk solution of colloids is dilute and only the translational 
entropy is relevant. Furthermore, the concentration of free chains and micelles of any size 
are also sufficiently dilute as to neglect interactions among them. 

Let c(m,p) be the number concentration of micelles with p arms and m colloids. Hence, 
the total free energy of the solution containing all kinds of aggregates is 



p 1//a 



VkT 



E ( c ( m 'P) ln + ( F ( m 'P) - /re/CO)) + c/m ^ - c f f s {T) 

m=Q p=l ' ' 

(1/a oo \ / 1/a oo \ 

c - E E mc ( m 'P) ~ c f) + ^ \ Cb ~ E ^2P c ( m 'P) ( 56 ) 
m=0 p=l I V m=0 p=l J 

where the last two terms fix the total amount of colloids and copolymer chains in the system, 
respectively. V is the volume of the system, c = N /V and = N b /V, where N and N b 
are, respectively, the total number of colloids and copolymers in the system. Furthermore, 
f s (T) stands for lnv s /A 3 s . In turn, f re f(T) = \nv/X 3 (see appendix IB1. Minimization of this 
free energy with respect to c(m,p) gives the equilibrium distribution of the aggregates by 
their size p and number of adhered colloids, m. That is, 
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vc(m,p) = MO, l)r(v s c f ) m ex V (- [F{m,p) - f ref - p(F(0, 1) - f ref ) - m(F(l, 0) - /,)]) 

(57) 

where the Lagrange multipliers /^o and Hb have been expressed through the concentrations of 
unimers of each species, c/ for free colloids and c(0, 1) for free and clean block copolymers, 
which can be both considered as the control parameters. 

To obtain the relationship between volume fractions of the species, we assume that there 
is no volume of mixing in the system. Then, the equilibrium distributions are given by 

<P(m,p) = P (0, 1)07 — x 

exp [- (F(m,p) - f ref -p(F(0, 1) - f ref ) - m(F(l,0) - /,))] (58) 

where <j)(m,p) is the volume fraction of a micelle of p arms with its adhered m stickers, while 
0(0, 1) and 4> s are, respectively, the volume fractions of free polymer and free colloid. 

To determine the effect of the presence of the adhered colloids on the micellization prop- 
erties, in the following we shall focus on the behavior of the function, 



U(m,p) = In 



c(m,p) 

= -(p- l)ln(c(0, l)v) -m\n(c f v s ) + F*{m,p) -pF* {0,1) (59) 

whose minima coincide with the maxima of c(m,p) and, therefore, gives us the average size 
of the micelle. Here F*(m,p) = F(m,p) — f re f for any p < 1 and we consider no energy 
associated with colloids. 

The convenience of this function is the possibility of naturally defining the critical ag- 
gregation concentration (cac), i.e., when the concentration of unimers and aggregates is of 
the same order, this potential is close to 0, while in the absence of the aggregation Q(m,p) 
goes to infinity. For convenience, we analyze Q(m,p) as a function of p for fixed m/p, to see 
how the presence of the adhered colloids change the properties of the "renormalized" block 
copolymer with m colloids per branch. In Figure [7] we make this analysis for different values 
of m/p. As we see the presence of adhered colloids induces the micellization and the growth 
of the micelles as m/p increased. 

We have focused our analysis on the regime where colloids are rather scarce along the 
polymer chains and the polymer concentration remains small in the associated micelles. 
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FIG. 7. Variation of 0(m,p) as a function of p for different values of m/p. cj has been set to e 10 
and c(0, 1), to e~ 30 , 7 = 0.7, L c = 50 and L = 1000. Parameters y = ^/kT^-du+i = 25 and 
a = A/L = 10 -4 . Notice that for c(0, 1) is below the CMC in the absence of colloids. The presence 
of colloids favors the micellization as well as the formation of larger aggregates as larger is m/p. 

Under these conditions the obtained micelle coronas are in the critical excluded volume 
regime and the bare interaction between effective monomers (sensitive to the presence of 
adsorbed colloids) is unimportant. Association is then favored by the effective shortening 
of the arms. However, at least close to the (small) core of the micelle, there may exist a 
region where the polymer is in the mean-field regime and therefore be sensitive to the bare 
magnitude of the excluded volume interaction. The effect of colloids on micellization then 
may be inverted there, that is, that colloid adhesion does not favor micellization. This regime 
deserves further consideration, although it lies beyond the scope of this work. Finally, if 
the chains are heavily loaded in colloids and lose flexibility the excluded volume interaction 
drops and micellization is then favored by the colloids. 



V. CONCLUSIONS 

A scaling theory is constructed aiming at the study of the conformations of hydrophilic 
polymers interacting with small colloidal particles that can reversibly stick onto a polymer 
backbone. Different geometries have been considered, namely linear polymer chains, star 
polymers, as well as coronas of block copolymer micelles. Unlike related analyses, we have 
consistently taken into account the integration of the position of the colloid over the chain 
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backbone, due to the reversible adhesion of the former onto the polymer. 

We find that small colloids adhered to a given polymer induce formation of independent 
loops of minimal size on each colloid, which can slide along the chain backbone, thus increas- 
ing the entropy gain in comparison with static loops of large size. This result is independent 
of the geometry of the polymer as well as of the nature of the interactions. The sliding 
degree of freedom of such small loops provides the system with an additional entropic con- 
tribution proportional to the length of the polymeric branch. More complex structures such 
as combined loops, ailed loops, or cactus-like structures are not favorable because of the 
lack of such a sliding degree of freedom. This contribution has been disregarded in previous 
analyses in similar systems that rely on the scaling of networks of fixed topology. This is 
one of the main results of this work. 

We have applied these concepts to the equilibrium between colloids and star polymers 
finding that the saturation of the polymer by the colloids is a crossover process, rather 
than a phase transition, except for the case of a star with infinitely long arms, where the 
crossover turns into higher order phase transition in the sense of Ehrenfest, that is, with a 
finite discontinuity in the third derivative of the free energy. Moreover, the analysis of the 
micellization problem in the presence of colloids agrees well with the observed experimental 
fact that the adhesion of colloids favors both, the micellization process, as well as the growth 
of larger micelles. 
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Appendix A: Ailed loop 

To illustrate this fact, we can for example calculate the partition function of ailed loops 
created by the presence of two stickers in the chain, as shown in Figure [3j Similar structures 
occur in knots formed by slip link contacts^. 

Denoting by a the length of the common section of the two ailed loops, whose total sizes 
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are n\ and n 2 respectively, and I is the size of the shortest tail, we can write the partition 
function as Z ~ a u nx 1 ri 2 i l y {N — I — n x — n 2 — a) z . 

Without loss of generality we choose n\ < n 2 , since the loops are symmetrical. In 
analogy with the previous case it can be shown that the configurations with the loops fixed 
at a free end, when the short tail satisfies I < n\ < n 2 are subdominant. Then, we consider 
n\ < n 2 < I- If the common section is small, a < ri\ < n 2 < I we can fix the exponents. 
If all segments are of the same order of magnitude and comparable to the total size of the 
chain, one has 

u + xi + x 2 + y + z = 2ai + 2cr 4 - 2du (Al) 
If the size of the common segment decreases to one single monomer, a — > 1, 

xi + x 2 + y + z = 2ai + a 6 - 2du (A2) 
If one loop is negligibly small, a < U\ — > 1, then 

x 2 + y + z = 2o\ + (T4 — du (A3) 
When both loops are vanishing, a < ri\ < n 2 — > 1, 

y + z = 2o x (A4) 

And, finally, small tail is vanishing, 

z = 2o"i (A5) 

From these equations all exponents can be calculated and the total partition function takes 
this form 



N/2 A pn 2 

dl / dn 2 n a 2 A - dv \ dnml 6 -^" x 



/ daa 2ai - ati (N — l — ri\—n 2 — a) 2 * 1 
Jo 



(A6) 



After integration, we find 
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Z ~ £2(<r<i-dv)+3 N 2ai 



(A7) 



The same scaling is obtained if one considers the other two possibilities, namely, n\ < n 2 < 
a < I and rii < a < n 2 < I and we have omitted its explicit calculation. As we have 
seen, the entropy favors the smaller size the possible for the loops. For this particular case 
this is obtained if both loops tend to shrink to its minimal possible size determined by the 
cut-off A. However since by construction loops are ailed, they are forced to move together 
along the chain backbone, therefore its entropy is a factor 1/N smaller than if they could 
independently move. This kind of reasoning can be extended to more complex topologies 
with the same result. 

If the sliding link is formed by small associating colloidal particles that can attach and 
detach from the chain at any point, the formation of individual loops moving separately 
along the chain is more favorable compared to such combined loops with common segments. 
Using the same arguments one can see that the partition function of individual loops ( [17]) 
dominates any other structures with the same number of stickers (combined loops, cactus 
structures of loops growing on top of each other, etc.). 

Appendix B: Calculation of the entropy of the micelles 

Self-assembly of polymers into finite size micelles is accompanied by the entropy changes 
due to the confinement of the polymers assembled into the micelle, which moves as a single 
object. Here, we estimate the entropic contribution due to this confinement, which involves 
the integration of the kinetic degrees of freedom, as well as the translational degree of 
freedom of the polymers as a whole. We recall that the conformation degrees of freedom 
have been explicitly considered in section [TT1 We are implicitly considering micellization as a 
two-state situation in which polymers are either free or aggregated into a micelle containing 
p polymers. 

First, we calculate the entropy change due to the association of the p polymers into a 
micelle in a system containing n of such micelles, to later introduce the effect of the presence 
of m stickers in each of them. Therefore, the partition function of n non-interacting micelles 
comprised of p undiscernable polymers, each of length N, is written as 
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1 r d'ft . . . (Tf Np yr f ~f\ 

p\J (2TihfP N ll eXP V 2m p kT 



d r i . . . d r p exp 



(Bl) 



kT J 

where m p is the mass of a monomer, U is the interaction potential between polymers forming 
the micelle (undefined). The integration is over all positions of the monomers r, and their 
momenta p. The micelles are distinguishable only by their composition. We recall that the 
integration over the translational degrees of freedom is constrained to the polymers to move 
together as an aggregate, thus the separation between any pair of monomers does not have 
to exeed the overall size of the aggregate. 

Equation (1B1I) can be rewritten in terms of the partition function of a single micelle of p 
polymers, Z p , according to Z = [Z p ] n , where Z p is written as 



Z p — — 



The momentum integral can be readily calculated to give a term inversely proportional 
to Ap pN , where A p is the de Broglie's length of the monomer (thus, independent of N) 
Ap = h 2 1 {2nm p kT) , and h, Planck's constant. In the spatial integration we will separate 
three contributions, namely, the overall translation of the aggregate, related to the center 
of mass motion, the relative motion of the polymers inside the aggregate, together with the 
displacement of the monomers of every polymer. The latter gives a contribution proportional 
to lp~^ per chain, where l p is a length of the order of the radius of gyration of the chain, 
which we will consider a molecular parameter. This contribution has effectively accounted 
for in the free energy term of the micelle that we have previously calculated (cf. sec. HB . 
Thus, due to the fact that in sec. [Til the entropy terms are given in a dimensionless form, 
the factor l p ^ simply provides the appropriate dimensionality of the expression. As far 
as the second contribution is concerned, the integration is limited to a distance between any 
pair of polymers of the order of the size of the aggregate, due to our initial hypothesis of 
the two-state system (free and confined), and the description that we are doing here of the 
precisely confined state. Therefore, this second contribution yields a term proportional to 
V^ 1 , where V agg stands for the volume of the aggregate. The remaining integration over the 
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position of the center of mass of the aggregate itself yields one additional factor proportional 
to the overal volume V of the system. Therefore, the final result reads 



Z p ~ 



yyp-lyip{N-\) 



p\Al pN 



(B3) 



Hence, we can write this entropic contribution to the free energy of the aggregate as 



F 
kT 



— In Z ~ n In 



n 

Ve' 



A* 



n(p — 1) In 



P 



-A 2 



— nln 



P 



(B4) 



where A = fl^~ x summarizes only molecular parameters. The first term is for the 
contribution of a ideal solution of micelles. The second and the third terms stand for the 
entropy penalty derived from keeping together p polymers into an aggregate. The third 
term is subdominant in the limit p > 1, If p — 1 this expression should recover the free 
energy of a dilute solution of independent polymers. Thus, the last term, in addition, 
should be regarded as an artifact of the distinction we have made between polymers inside 
the aggregate and outside the aggregate. By consistency, the last term should be dropped 
to match that mentioned limit, provided that in the more interesting situation p ^> 1 is 
subdominant. Notice that the calculations done in sec. HTItake into account the entropy due 
to the conformational changes of a system with p arms, while here we are calculating the 
entropy penalty of putting these p arms together into an aggregate, understanding that the 
aggregation process is a physical process. Therefore, V agg should be regarded as the volume 
of the core of the micelle and not that of the whole aggregate, including the corona. 

Similar arguments have been used in Refs.— ^ to derive the entropy penalty for the 
formation of the micelle, obtaining an expression similar to the second term in our eq. (1B4I) . 
The main difference comes from the presence in our expression of the A 3 term. Therefore, 



the entropic contribution of Ref. 
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should be regarded as with respect to a reference state. 
Effectively, let us introduce a volume Vq, where Vq = Nv, v being the volume of the monomer, 
into the logarithmic term, and substract it afterwards. The free energy can be written as 



kT 



nln 



n 

leV 



Vn 



+ n(p — 1) In 



P 



-V 



a-gg 



— np 



ref 



kT 



where f re f is the free energy per polymer in a reference state, defined as 



fref 
kT 



hi 



Vo 
A 3 



(B5) 



(B6) 
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Eq. (1B5I) is exactly the same expression as in Ref. 



37 



In the case of the colloids, the situation is different, for the entropy of the adhered colloid 
has already taken into account in the derivation of Z(m,p). Notice that this expression 
accounts for the energetic contribution exp(— em), the translational entropy due to the 
colloid motion along the chain backbone, as well as the indiscernibility of the colloids. 
Therefore, we are only facing the contribution due to the chain aggregation. 

Finally, for convenience, we introduce the volume of the monomer v instead of Vq into 
the expression (IB4I) . Then, eq. IB4I can be rewritten as 



7^ =-lnZ 
kT 



nln 



n 

V? 



+ n(p — 1) In 



V 



npfref(T) 



(B7) 



where f re f(T) = lni;/A 3 which is only dependent on the temperature and molecular param- 
eters, and, in addition, is independent of the size of the polymer. We use this function in 
the construction of the free energy in eq. ( |56|) for the systems of micelles of different size 
and number of colloids in equilibrium. 



Therefore, from this expression we can identify that the contribution that yields F con f is 
given by the second term in IB7[ that is 



nF conf /kT = n(p - 1) In 



P 



[VaggZ 



(B8) 



As it has been mentioned, the other contributions are explicitly accounted for in the 
proposed expression for the free energy, eq. fl56|) . 
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